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Precise measurements of inclusive multiplicities of charged pions and kaons produced in e + e~ 
annihilation have been carried out at a center-of-mass energy of \fs = 10.52 GeV. The measurements 
were performed with the Belle detector at the KEKB e + e~ collider using a data sample containing 
113 x 10 6 e + e~ — > qq events, where q — {u,d, s,c}. We present charge-integrated multiplicities 
1/iVevts x dN h ±/dz for h ± = { 7 r ± ,if ± } as a function of the relative hadron energy z = 2E h /^/s 
from 0.2 to 0.98 (0.97) for ti ± (K ). The combined statistical and systematic uncertainties for 
tt* (A' ± ) are 4% (5%) at z ~ 0.6 and 18% (51%) at z ~ 0.9. The multiplicities are the first 
measurements of the z-dependence of pion and kaon multiplicities for z > 0.7 as well as the first 
precision multiplicity measurements at a center-of-mass energy far below the Z° resonance used by 
the experiments at LEP and SLC. 

PACS numbers: 13.66.Bc 



The fragmentation function (FF) D^- g (z,Q 2 ) 
parametrizes the transition of quarks q, antiquarks q and 
gluons g into the color-neutral hadron h in high-energy 
particle collisions as a function of the four-momentum 
transfer squared, Q 2 , in the collision and the relative 
energy carried by the final state hadron, z = 2Ef l /^/s. 
FFs are process-independent and are extracted through 
global analysis of hadron multiplicity measurements 
obtained from inclusive hadron production in e + e~ and 
pp reactions and from semi-inclusive hadron production 
in IN reactions at various center-of-mass energies (see 
Refs. [l|-[3j])- Dokshitzcr-Gribov-Lipatov-Altarelli-Parisi 
(DGLAP) evolution equations |ij derived from Quantum 
Chromodynamics (QCD) are used to jointly analyze data 
sets taken at different Q 2 . The Q 2 evolution mixes quark 
and gluon degrees of freedom and enables the extraction 
of quark and gluon FFs through the analysis of precise 



hadron multiplicity data sets. In e + e~ annihilation, the 
availability of precision data sets was previously limited 
mostly to results from experiments (e.g., Refs. 0-0]) at 
LEP and SLC at the energy scale of the Z°-boson mass. 
The large data samples available at Belle complement 
existing results on hadron multiplicities with precise 
measurements at Q 2 = (10.52 GeV) 2 . The multiplicities 
reported here are expected to improve the precision of 
FFs and the QCD-based determination of the nucleon 
quark- and gluon-structure from semi-inclusive deeply 
inelastic measurements at CERN, DESY and Jefferson 
Laboratory and inclusive proton-proton measurements 
at CERN and BNL. 

This analysis is based on a 68.0 fb^ 1 data sample col- 
lected with the Belle detector at the KEKB asymmetric- 
energy e+e- (3.5 on 8 GeV) collider @ operating 60 McV 
below the T(45) resonance at = 10.52 GeV. The 
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Belle detector is a largc-solid-angle magnetic- spectrom- 
eter that consists of a silicon vertex detector (SVD), a 
50-layer central drift chamber (CDC) and an array of 
aerogel threshold Cherenkov counters (ACC). The ACC 
detector is surrounded by a barrel-like arrangement of 
time-of-flight scintillation counters (ToF) and an electro- 
magnetic calorimeter (ECL) comprised of CsI(Tl) crys- 
tals, located inside a superconducting solenoid coil that 
provides a 1.5 T magnetic field. The iron flux-return yoke 
surrounding the coil is instrumented to detect K\ mesons 
and to identify muons (KLM). The detector is described 
in detail elsewhere 0. The Belle Monte Carlo (MC) 
simulations are performed with the PYTHIA 6.205 [llj 
event generator and a GEANT3 [ijj detector simulation. 

An event must have at least three charged tracks and 
a heavy-jet mass (the sum of all reconstructed parti- 
cle masses in the jet having the largest reconstructed 
mass), above 1.8 GeV/c 2 or above 25% of the recon- 
structed visible energy. A reconstructed visible energy 
of more than 7 GeV is required to reduce contamina- 
tion from t + t~ events. From these events, we retain 
tracks with a laboratory-frame momentum pi a b of more 
than 500 MeV/c, a scattering angle 9\ a h (relative to 
the beam axis) within the central detector acceptance 
(-0.511 < cos6»i ab < 0.842), at least three SVD hits, and 
a point of closest approach to the e + e~ interaction point 
of under 4 cm along the beam axis and under 1.3 cm in 
the perpendicular plane. 

Using likelihood ratios based on dE/dx measurements 
of the CDC, signals from the ACC, ToF measurements, 
ECL cluster energies as well as KLM signals, charged 
particles are identified as pions, kaons, protons, electrons 
or muons. The initially-measured track yield Nj iineas for 
particle type j is binned in z with width 0.01, starting 
from 0.2; negatively- and positively-charged tracks are 
treated separately. In order to obtain final multiplici- 
ties, several corrections are applied, as described below. 
Statistical uncertainties are propagated through all cor- 
rections and systematic uncertainties from all corrections 
are added in quadrature. 

The true yield iV^true for particle type i is obtained 
from the measured yields via 

N^tmc = M U ljV J^cas, (1) 

j 

where A/ -1 is the inverse of the 5x5 particle iden- 
tification (PID) probability matrix M, whose diago- 
nal and off-diagonal elements represent PID efficien- 
cies and misidentification probabilities, respectively. As 
the elements of M depend on the laboratory-frame 
momentum and polar angle, an M matrix is con- 
structed in each of the 17 x 9 bins with p\ a ^ bound- 
aries [0.5,0.65,0.8, 1.0, ...,3.0,3.5,4.0,5.0,8.0) GeV/c 
and cost9 lab boundaries [-0.511,-0.300,-0.152,0.017, 
0.209, 0.355, 0.435, 0.542, 0.692, 0.842). These values were 
chosen to obtain similar track yields in all bins. The 



elements of M, called "efficiencies" henceforth, are ob- 
tained mostly from experimental data. For pions and 
kaons, the decay D*+ -+ tt+ ow + D° -> ^+ ow + (K'tt+J 
and its charge-conjugate are used. Similarly, proton 
PID efficiencies are extracted from A — > p + + tt~ and 
charge-conjugate decays and electron and muon efficien- 
cies from Jji\) — > e + e~, /i + /J~ decays. The procedures 
for calculating PID efficiencies from decays are rather 
similar and are described here for the pions in the D* 
channel. The reconstructed invariant mass distribution 
TO 7r B iow(7Tfa 5 t^") — TO (7Tf aB t-fO f° r candidate triplets of tracks is 
fitted with a threshold background function and a mod- 
ified Crystal Ball [IH function for the signal. The initial 
yield of true pions in the sample can be obtained from the 
fitted yield of D*'s in the invariant mass distribution be- 
fore applying any PID selection criteria. Performing the 
fit on the same track sample but with additional PID se- 
lection criteria for pions, kaons, protons and leptons, and 
dividing by the initial yield of true pions, gives pion PID 
efficiencies. All fit uncertainties are propagated through 
the efficiencies, accounting for correlations. When no sig- 
nal can be fitted, the efficiency is set to zero using the 
largest possible signal size as the uncertainty. Several 
matrix elements, especially the ones for leptons due to 
the high J/ip mass relative to the quark energy, cannot 
be filled only from data-driven decay studies. Missing ef- 
ficiencies are estimated once by extrapolating the trend 
of extracted probabilities with MC information and a sec- 
ond time by utilizing MC probabilities only, which pro- 
duces two alternative sets of PID matrices; both are used. 
Examples of pion efficiencies are shown in Fig. [TJ Gen- 
erally, the diagonal matrix elements have values around 
90% with the exception of kaon and proton efficiencies 
that drop at the highest momenta to about 70% and 
50%, respectively. The most prominent misidentification 
probabilities are it — > K (up to 15%), K — > 7r (up to 
20%) and p — > K (up to 50%). For the correction ac- 
cording to EqJTJ both sets of PID matrices are inverted 
using singular value decomposition tools [l4[ . Uncertain- 
ties on the elements of A/ -1 are obtained by statisti- 
cal analysis of large samples of inverted PID matrices 
whose elements are varied before inversion according to 
their uncertainties. Finally, these uncertainties are prop- 
agated as systematic uncertainties to the PID-corrected 
track yields. Measured yields iVj jmeas are corrected with 
each set of alternative matrices. The final PID-corrected 
yields are calculated from the average of the two correc- 
tions. The difference between the corrections is assigned 
as additional systematic uncertainty. The PID correc- 
tion changes pion (kaon) yields by around +10% (—20%) 
at low z and by around +5% (+10%) at high z. The 
dependence on z of the changes is due to the kinematic 
dependence of the extracted PID efficiencies. 

After the basic hadronic event selection, experimen- 
tal data still contain hadron contributions from events 
other than e + e~ — > qq. A MC sample containing all non- 
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FIG. 1. (color online). Samples of extracted pion PID ef- 
ficiencies (green triangles) and pion misidentification proba- 
bilities to kaons (blue squares), protons (red circles), muons 
(purple crosses) and electrons (green diamonds) for pions with 
0.4 < cos^iab < 0.5 as a function of the pion laboratory-frame 
momentum pi a b- Electron probabilities are small and overlap 
with muon probability symbols in the plot. 



qq physically possible event topologies in e + e annihila- 
tion at Q 2 = (10.52 GeV) 2 is assembled from dedicated 
QED (BHLUMI [3 generator for Bhabha-scattering, 
KKMC generator for [i and t pair production) and 
two-photon MC data. Hadron yields from such non-gg 
events passing the hadronic event selection are scaled 
to the luminosity of the measurement sample and sub- 
tracted from the PID-corrected yields. Contributions 
from events other than e + e~ — »• qq are found not to be 
significant at low z. At the maximum z, about 30% of the 
reconstructed pions and 10% of the reconstructed kaons 
are created in non-QCD rr and two-photon events. Sta- 
tistical uncertainties are propagated as systematic uncer- 
tainties on the corrected yields. 

The reconstructed momenta of charged tracks and 
hence the reconstructed values of z are smeared around 
the true values because of the momentum resolution. 
Similar to the particle misidentification correction, this 
z-bin migration can be described by a matrix. While a 
z cut-off of 0.2 is applied to the final multiplicities, the 
smearing matrix needs to be evaluated for lower z to take 
into account smearing into or out of the selection region. 
Smearing matrices are extracted down to z = 0.08 for 
pions and z = 0.12 for kaons from the MC and detector 
simulation. The matrices are inverted and uncertainties 
are propagated with the same techniques as for the PID 
correction. Due to limitations of the method, the correc- 
tion is only applied to yields below z = 0.98 for pions 
and z = 0.97 for kaons and so imposes an upper limit on 
z for the final multiplicities. All uncertainties on inverse 
matrix elements are propagated to the corrected yields 
as systematic uncertainties and dominate the total un- 
certainties for high z. The smearing correction causes no 



significant change in the measured yields because of the 
good momentum resolution, small bin width and mostly 
symmetric smearing to lower and higher values of z. The 
momentum resolution of the detector simulation has been 
compared to the resolution extracted from experimental 
cosmic ray data. The small differences arc assigned as an 
additional systematic uncertainty. 

Some hadrons are lost or gained through decays in 
flight. In addition, hadrons are lost due to interaction 
with detector material while others are created in such 
interactions. Finally, the track finding and fitting algo- 
rithm may find fake tracks or fail to reconstruct true 
tracks because of detector inefficiencies or occupancy ef- 
fects. Limitations in the MC do not allow us to distin- 
guish these processes. Thus, a joint correction is per- 
formed for all processes in each z bin, extracting ra- 
tios between reconstructed and generated particle yields. 
Concerning decay- in-flight (DIF), the QCD framework 
used by global analyses usually does not account for 
hadrons created in weak decays, which are nevertheless 
present in measured multiplicities. Correction factors to 
remove pions and kaons produced in weak decays (e.g., 
from t leptons, kaons and D mesons, A and heavier 
baryons) would be purely MC-dcpcndcnt. To be consis- 
tent, all weakly produced pions and kaons are included 
in the results, and all decayed pions and kaons are recov- 
ered. For the final multiplicities, we provide supplemen- 
tary fractions of pions and kaons originating from strong 
and weak decays [15|, extracted from our MC data. At 
low to medium z, the fraction of weakly-produced pions 
(kaons) is 30% (50%) but vanishes towards high z due 
to phase-space limitations. The implementation of DIF 
in the MC is tested by analytical calculations and found 
to be consistent. Correlation effects between the PID 
and DIF corrections are studied; the observed correla- 
tions are applied as additional systematic uncertainties. 
The effect of hadronic interaction modeling on particle 
yields is estimated by comparing efficiencies from the de- 
fault FLUKA with the GEISHA ^ package and 
with no hadronic interactions. The small remaining dif- 
ferences are assigned as systematic uncertainties. Fake- 
and multiple-track reconstruction is found to be negli- 
gible. Track reconstruction inefficiencies are at the few- 
percent level and are corrected for. The most significant 
track quality cut is the requirement for three or more 
hits in the SVD, which also introduces most of the high- 
z dependence of the efficiencies in this correction. The 
overall efficiencies remain above 85% for all z. Statis- 
tical and systematic uncertainties on the efficiencies are 
propagated as systematic uncertainties on the corrected 
yields. 

To obtain the best-measured reconstructed tracks, the 
analysis is limited to the barrel region of the detector. A 
correction is applied to recover the 4-7T physical particle 
yield. Corresponding efficiencies from the MC decrease 
from about 65% at low z and level out around 60% at high 
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z. The high-z behavior is consistent with the assumption 
that high- z hadrons follow the known 1 + cos 2 6 cms distri- 
bution of fragmenting quarks, where 9 cms is the scattering 
angle in the center-of-mass frame. Towards lower z, addi- 
tional transverse hadron momentum causes the cos 9 cms 
spectra to flatten and efficiencies to increase. The MC 
description of the cos 9 cms dependence in the fragmenta- 
tion process is tested with experimental data yields in a 
two-dimensional z and cos0 cms binning. Consistency is 
found within statistical uncertainties and no additional 
systematic uncertainty is assigned. Statistical uncertain- 
ties on the efficiencies are propagated as systematic un- 
certainties on the corrected yields. 

Hard initial-state and final-state photon radiation 
(1SR/FSR) processes reduce the fragmentation energy 
scale f° r both (ISR) or for one (FSR) final-state 

quark. Therefore, experimentally measured yields con- 
tain a variety of fragmentation scales y/s/2. A theoret- 
ical deconvolution of ISR/FSR from measured yields is 
beyond the scope of this measurement. Instead, we keep 
the energy scales in the measurement sample within 0.5% 
of the nominal y/s/2 in order to remain below the scale 
resolution of the state-of-the-art next-to-leading order 
DGLAP evolutions. Corresponding fractions of hadrons 
from events with summed ISR/FSR photon energies of 
less than 0.5% x y/s/2 are extracted from the MC for 
each z bin. These fractions rise from 65% at low z to 
almost 100% at maximum z due to phase-space limita- 
tions for ISR/FSR. The fractions are applied bin- by- bin 
to the measured yields to exclude particles from events 
with large ISR/FSR contributions. Systematic uncer- 
tainties are assigned from the dependence of the hadron 
fractions on the chosen PYTHIA MC parameter sets and 
from accounting for the remaining scale variance in the 
sample. 

Efficiencies for all applied event selections are ex- 
tracted from MC simulations. Efficiencies for the 
hadronic event selection are close to unity at low z. They 
drop rapidly at larger z since both the track multiplic- 
ity and heavy jet mass requirements disfavor events with 
high-z tracks. In such events, the phase space for addi- 
tional particles and substantial remaining transverse mo- 
mentum is significantly reduced. The efficiencies drop to 
about 60% at z ~ 0.9 and to less than 10% for z ~ 1.0. 
Conversely, efficiencies for the visible energy requirement 
are above 95% at high z, but drop to about 75% for low- 
z particles. This can be understood from the increased 
likelihood for events with high-z tracks to deposit suf- 
ficient energy in the barrel part of the detector. The 
measured yields are corrected accordingly and statistical 
uncertainties are propagated as systematic uncertainties 
on the corrected yields. The distributions of event shape 
variables used in the event selection are compared be- 
tween MC and experimental data. Differences in these 
distributions causing variations in event selection efficien- 
cies are assigned as systematic uncertainties. 



After all corrections, the measured yields arc nor- 
malized to the total number of hadronic events in the 
measurement sample. The latter is calculated from the 
product of the total luminosity of 68.0 fb" 1 with the 
corresponding total hadronic cross-section. This cross- 
section is extracted from a measurement of the R frac- 
tion, R = a(e + e~ — > h + X)/<7q [2(|, and the calculated 
point-like cross-section <jq = 47rag M /(3s). Taking into 
account that about 35.4% of all events are excluded from 
the analysis due to large ISR/FSR contributions, a total 
normalization of (122.8 ± 3.0) x 10 6 e + e~ —J- qq events is 
obtained, where q = {u, d, s, c}. 

The resulting final charge-integrated pion and kaon 
multiplicities 1/N cvts x dN h ±/dz for ft ± = {n ± ,K ± } 
are displayed in Fig. [5J All multiplicity values and un- 
certainties are available in Ref. [15| . Combined relative 
statistical and systematic uncertainties remain below 6% 
up to fractional hadron energies z ~ 0.6, then rise to 
about 18% (pions) and 51% (kaons) at z ~ 0.9 and reach 
70% (pions) and more than 100% (kaons) at maximum 
z. As a test of all applied corrections, pion and kaon 
charge ratios Ri = N i -/N i + are fitted with a constant. 
The results, = 0.996 ±0.010 and R K = 1.000 ±0.011 
(with a combined statistical and systematic uncertainty), 
are consistent with 1.0 within the extracted uncertainties 
and indicate consistency of the performed systematic cor- 
rections. The precision of the present measurement is 





^k* (Statistical, Systematic Uncertainty) 
-4JK 1 (Statistical, Systematic Uncertainty) 
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FIG. 2. (color online). Final charge-integrated pion (green 
solid boxes) and kaon (blue empty boxes) multiplicities as a 
function of the fractional hadron energy z. Statistical uncer- 
tainties are shown by error bars, systematic uncertainties by 
error bands. 

systematics- limited for all z. Figure [3] shows different 
contributions to the systematic uncertainties for the fi- 
nal charge- integrated pion multiplicities in a stacked plot. 
At low z, the dominant systematic uncertainties arise 
from the initial/final state radiation correction. At high 
z, the momentum smearing, particle identification and 
DIF/rcconstruction corrections represent the largest con- 
tributions to the systematic uncertainties. An additional 
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constant overall normalization uncertainty of 2.4% is not 
shown. 
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FIG. 3. (color online). Stacked plot of the increase of rel- 
ative systematic uncertainties due to each correction step. 
The individual contributions are particle identification (solid 
gray), sample impurity (vertical lines), momentum smearing 
(oblique lines), decay- in-flight and reconstruction (horizontal 
lines), event /track selection (star hatch), acceptance (cross 
hatch) and initial/final state photon radiation (brick hatch). 
The relative statistical uncertainty is shown in red dashes. An 
additional normalization uncertainty of 2.4% is not shown. 



scribed above is lower or comparable to the precision 
achieved in the results reported here. In particular, sig- 
nificantly better precision than previous measurements 
at low energy scales is reached. 



"6 

■D 



10 7 
10 6 
10 5 
10" 
10 3 
10 2 



a (Sel.) for e *e -> ji=+X Multiplicities 





















OoooftpSlGeVf, 


10 ,! ) 








: ^ 






















10") ° 








* rv 




> ' 








■-» (J 


IGeV 






■ 










(»7xiS') 












to 6 ; 

• 












o GeV WOOD) ■»» 


7 0Q0Q } 










I0G e 


*<X1S0) 










^ - . "onan m „. 








■ • 




i i "*• , , , 







0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



FIG. 4. (color online). Charge-averaged pion multiplicities 
from selected measurements. For better visibility, all data sets 
are scaled by arbitrary numbers (indicated between parenthe- 
sesjf|- Statistical and systematic uncertainties are added in 
quadrature. 



The final pion and kaon multiplicities are compared 
to multiplicities from the Belle MC and other MC sam- 
ples generated with the PYTHIA 6.205 generator for 
e + e~ — > qq events, where q = {u,d,s,c}, at our energies. 
PYTHIA parameter sets as used by higher energy exper- 
iments such as CDF, ATLAS and ALEPH are used. The 
{u, d, s, c}-quark Belle MC is optimized together with the 
6-quark MC to describe low-z particles created in T(4S*) 
resonance production. Comparing to the Belle measure- 
ment, the {u, d, s, c}-quark MC significantly underesti- 
mates the data for z > 0.45. Multiplicities generated 
with the ALEPH parameter set overshoot Belle experi- 
mental data for all z. Multiplicities from CDF and AT- 
LAS parameter sets describe the Belle results within 20%. 
Figures of ratios of MC over experimental data multi- 
plicities for pions and kaons with statistical and system- 
atic uncertainties combined and propagated are given in 
Ref. [H. 

The final multiplicities are compared with identified 
multiplicities from LEP and SLC experiments as well as 
from other, lower energy e + e~ experiments. As an il- 
lustration, the pion multiplicities are compared to those 
from Refs. @, i, H in Fig. g] The correspond- 

ing plot for kaon multiplicities can be found in Ref. . 
The resolution in z is significantly improved for all com- 
pared multiplicities over most of the z range of this mea- 
surement. In addition, no other previous measurement 
probes the z dependence of multiplicities at high z > 0.7. 
The relative precision of the previous measurements de- 



Presentation style based on plots created by O. Bicbcl et al. 
shown in Ref. '27il 



In conclusion, multiplicities of identified charged pi- 
ons and kaons are measured over a broad range in z with 
0.2 < z < 0.98 (tt), 0.97 (K) with high relative precision. 
The analyzed data sample of 68.0 fb _1 has been accu- 
mulated at a center-of-mass energy of y/s = 10.52 GeV, 
rendering this measurement the first precision measure- 
ment far from LEP/SLC center-of-mass energies. The 
high statistics and good control of systematic uncertain- 
ties will give constraints on the z-dependence of hadron 
FFs at large z > 0.7 for the first time. 
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Supplemental Document to Manuscript 
Precision measurement of charged pion and kaon multiplicities in e + e~ annihilation at 

y/s = 10.52 GeV 



All multiplicity values and uncertainties, together with the fractions of hadrons from strong and weak decays are 
not confirmed as final at this time because of the ongoing manuscript review at PRL. The data will be confirmed as 
final after the manuscript has been accepted for publication by PRL. 

CHARGE-INTEGRATED PION MULTIPLICITIES 

Low z = 2Eh/\fs bin limit (inclusive, constant bin width Az = 0.01), multiplicity value l/N ev t B x dN^±/dz, 
statistical uncertainty, positive systematic uncertainty, negative systematic uncertainty. 
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CHARGE-INTEGRATED KAON MULTIPLICITIES 

Low z = 2Eh/\fs bin limit (inclusive, constant bin width Az = 0.01), multiplicity value 1/N ev t s x dN^± /dz, 
statistical uncertainty, positive systematic uncertainty, negative systematic uncertainty. 
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FRACTIONS OF PIONS PRODUCED IN STRONG AND WEAK DECAYS 

Low z = 2Eh/\fs bin limit (inclusive, width Az = 0.01), fraction from strong decays, uncertainty of fraction from 
strong decays (stat.), fraction from weak decays, uncertainty of fraction from weak decays (stat.). 
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FRACTIONS OF KAONS PRODUCED IN STRONG AND WEAK DECAYS 

Low z = 2E] 1 /t/s bin limit (inclusive, width Az = 0.01), fraction from strong decays, uncertainty of fraction from 
strong decays (stat.), fraction from weak decays, uncertainty of fraction from weak decays (stat.). 
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COMPARISON OF MEASUREMENT RESULTS WITH MONTE CARLO MULTIPLICITIES 

The final charge-integrated pion and kaon multiplicities, 1/N evts x dN h ± /dz with z = 2Eh/yfs and = {n^ , K^}, 
are compared with multiplicities from the Belle MC and other MC samples generated with the PYTHIA 6.205 
generator for e + e~ — > qq events, where q = {u,d, s, c}, at Belle center-of-mass energy. PYTHIA parameter sets 
as used by high-energy experiments like CDF, ATLAS and ALEPH are studied. Figure [T] shows ratios of MC 
over experimental data multiplicities for pions (left) and kaons (right) with statistical and systematic uncertainties 
combined and propagated. 




FIG. 1: Ratios between PYTHIA Monte Carlo and experimental data multiplicities for pions (left) and kaons (right), for the 
Belle MC (circles), ALEPH (boxes) and CDF/ ATLAS parameter sets (triangles). 



COMPARISON OF KAON MULTIPLICITIES WITH WORLD MULTIPLICITY DATA 



The final charge-integrated kaon multiplicities 1/N cvts x dN K ±/dz with z = 2Eh/^/s are compared with kaon 
multiplicities from LEP and SLC experiments as well as from other, lower energy e + e~ experiments from Refs. 0- 
in Fig. |3 
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FIG. 2: Charge-averaged kaon multiplicities from selected measurements. For better visibility all data sets are scaled by 
arbitrary numbers (indicated between parentheses). Statistical and systematic uncertainties are added in quadrature. 
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